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Abstract. We performed an extended analysis of the parameter space for the interaction of the Magellanic System with the 
Milky Way (MW). The varied parameters cover the phase space parameters, the masses, the structure, and the orientation of 
both Magellanic Clouds, as well as the flattening of the dark matter halo of the MW. The analysis was done by a specially 
adopted optimization code searching for the best match between numerical models and the detailed H I map of the Magellanic 
System by Briins et al. 1 2005 I. The applied search algorithm is a genetic algorithm combined with a code based on the fast, but 
approximative restricted N-body method. By this, we were able to analyze more than 10* models, which makes this study one of 
the most extended ones for the Magellanic System. Here we focus on the flattening q of the axially symmetric MW dark matter 
halo potential, that is studied within the range 0.74 < q < 1.20. We show that creation of a trailing tail (Magellanic Stream) and 
a leading stream (Leading Arm) is quite a common feature of the Magellanic System-MW interaction, and such structures were 
modeled across the entire range of halo flattening values. However, important differences exist between the models, concerning 
density distribution and kinematics of H I, and also the dynamical evolution of the Magellanic System. Detailed analysis of 
the overall agreement between modeled and observed distribution of neutral hydrogen shows that the models assuming an 
oblate (q < 1.0) dark matter halo of the Galaxy allow for better satisfaction of HI observations than models with other halo 
configurations. 
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1. Introduction 

The idea of dark matter (DM) was introduced by 
Zwicky 0933> . His dynamical measurements of the mass-to- 
light ratio of the Coma cluster gave larger values than those 
known from luminous parts of nearby spirals. That discrepancy 
was explained by the presence of DM. Ostriker et al. (I1974> 
proposed that DM is concentrated in a form of extended 
galactic halos. Analysis of rotation curves of spiral galaxies 
(Bosma 119811 Rubin & Burstein I1985> denotes that their 
profiles cannot be explained without presence of non-radiating 
DM. Hot X-ray emitting halos have been used to estimate total 
galactic masses (McLaughlin 1199^ . Corresponding mass-to- 
light ratios exceed the maximum values for stellar populations, 
and DM explains the missing matter naturally. The presence 
of DM halos is expected by the standard CDM cosmology 
model of hierarchical galaxy formation. The classical CDM 
halo profile (NFW) is simplified to be spherical. However, 
most CDM models expect even significant deviations from the 
spherical symmetry of DM distribution in halos. The model of 
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formation of DM halos in the universe dominated by CDM by 
Frenk et al. ( I1988t produced triaxial halos with a preference 
for prolate configurations. Numerical simulations of DM halo 
formation by Dubinski & Carlberg (1991) are consistent with 
halos that are triaxial and flat, with (c/a) = 0.50 and (b/a) - 
0.71. There are roughly equal numbers of dark halos with 
oblate and prolate forms. 

Observationally, the measurement of the shape of a DM 
halo is a difficult task. A large number of varying techniques 
found notably different values, and it is even not clear if the 
halo is prolate or oblate. Oiling & Merrifield (20001 use two 
approaches to investigate the DM halo shape of the Milky Way 
(MW), a rotation curve analysis and the radial dependence of 
the thickness of the H I layer. Both methods lead consistently 
to flattened oblate halos. 

Recently, the nearly planar distribution of the observed 
MW satellites, which is almost orthogonal to the Galactic 
plane, raised the question if they are in agreement with cos- 
mological CDM models (Kroupa et al. 2005 1 or if other ori- 
gins have to be invoked. Zentner et al. i2005j claim that the 
disk-Uke distribution of the MW satelUtes can be explained, 
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provided the halo of the MW is sufficiently prolate in agree- 
ment with their CDM simulations. On the other hand, it is 
not clear if there exists a unique prediction of the axis ratios 
from CDM simulations, as the scatter in axis ratios demon- 
strates (Dubinski & Carlberg 1991). Based on ACDM simu- 
lations, Kazantzidis et al. (2004i emphasize that gas cool- 
ing strongly affects halo shapes with the tendency to produce 
rounder halos. 

Another promising method to determine the Galactic halo 
shape are stellar streams because they are coherent structures 
covering large areas in space. Thus, their shape and kinematics 
should be strongly influenced by the overall properties of the 
underlying potential. A good candidate for such an analysis is 
the stellar stream associated with the Sagittarius dwarf galaxy. 
By comparison with simulations, Ibata et al. (2001 1 found that 
the DM halo is almost spherical in the galactocentric distance 
range from 16 to 60 kpc. Helmi (2004i studied the Sagittarius 
dwarf stream and argues for a non-spherical shape of the 
Galactic halo. However, Helmi (2004i also warned that the 
Sagittarius stream might be dynamically too young to allow 
for constraints on the halo shape. Johnston et al. (20051 argue 
that the analysis of the orbital planes of the Sagittarius dwarf 
galaxy prefers oblate over prolate models. In flattened oblate 
halos, the non-polar satellite orbits tend to become co-planar 
and contribute to the disk formation. In the case of Monoceros 
tidal stream discussed by Penarrubia et al. J2005t . the best- 
fitting models predict a significantly oblate halo. 

In this paper we use the Magellanic Stream to derive 
constraints on the halo shape of the MW. The basis are the 
new detailed HI observations of the Magellanic System (in- 
cluding the Large Magellanic Cloud (LMC) and the Small 
Mageflanic Cloud (SMC)) by Bruns et al. i2005> . As rem- 
nants of the LMC-SMC-MW interaction, extended struc- 
tures connected to the System are observed. Among them, 
the Magellanic Stream - an HI tail originating in be- 
tween the Clouds and spreading over x 100° of the plane 
of sky - has been a subject of investigation for previ- 
ous studies (see, e.g., Fujimoto & Sofue 1976; Lin&Lynden- 
Bell [T9771 Murai& Fujimoto 1980 Heller & Rohlfs l^ 
Gard iner & Noguchi [1996 ; Bekki&Chiba 2005; Mastropietro 
et aI. I2005t . Due to the extended parameter space related to the 
interaction of three galaxies and also due to the high compu- 
tational costs of fully self-consistent simulations, simplifying 
assumptions were unavoidable. Many simulations neglected 
the self-gravity of the individual stellar systems by applying 
a restricted N-body method similar to the method introduced 
by Toomre & Toomre (1972i. None of these simulations con- 
sidered the self-gravity of all three galaxies. Often only one 
galaxy is simulated, including its self-gravity by means of a 
live disk and halo, whereas the other two galaxies are taken into 
account by rigid potentials of high internal symmetry. That is, 
none of the simulations so far adopted a live dark matter halo 
of the MW, but they applied (semi-)analytical descriptions for 
the dynamical friction between the Magellanic Clouds and the 
MW. Also, a possible flattening of the MW halo has not been 
considered. Having the numerical difficulties in mind, it is not 
surprising that a thorough investigation of the complete param- 
eter space was impossible. 



Modeling observed interacting galaxies means dealing with 
an extended high-dimensional space of initial conditions and 
parameters of the interaction. Wahde and Theis ('1 999> 

introduced a genetic algorithm (GA) as a robust search method 
to constrain models of observed interacting galaxies. The GA 
optimization scheme selects models according to their ability 
to match observations. Inspired by their results, we employed a 
simple fast numerical model of the Magellanic System com- 
bined with an implementation of a GA to perform the first 
very extended search of the parameter space for the interac- 
tion between LMC, SMC, and MW. Here we present our re- 
sults about the MW DM halo flattening values compatible with 
most detailed currently available H I Magellanic survey (Briins 
et al. l2005l see Figs.^ andEJi. 

2. Magellanic Clouds and MW interaction 

2.1. Observations of the Magellanic System 

The MW, together with its close dwarf companions the LMC 
and SMC, forms an interacting system. Hindman et al. J1963t 
observed the HI Magellanic Bridge (MB) connecting the 
Clouds. Another significant argument for the LMC-SMC- 
MW interaction was brought by Wannier & Wrixon ( 1 972J and 
Wannier et al. J1972t . Their H I observations of the Magellanic 
Clouds discovered large filamentary structures projected on the 
plane of the sky close to the Clouds, and extended to both 
high negative and positive radial LSR velocities. Mathewson 
et al. (1974 1 detected another HI structure and identified a 
narrow tail emanating from the space between the LMC and 
SMC, spread over the South Galactic Pole. The tail was named 
the Magellanic Stream. A similar HI structure called the 
Leading Arm extends to the north of the Clouds, crossing the 
Galactic plane. A high-resolution, spatially complete HI sur- 
vey of the entire Magellanic System done by Briins et al. (120051 
see Figs.n and|2} gives detailed kinematic information about 
the Clouds and the connected extended structures. It indicates 
that the observed features consist of the matter torn off" the 
Magellanic Clouds and spread out due to the interaction be- 
tween the LMC, SMC, and MW. 

2.2. Modeling of the Magellanic System 

Toorrtre & Toomre (1972i have shown the applicability of re- 
stricted N-body models on interacting galaxies. In restricted 
N-body simulations, gravitating particles are replaced by test- 
particles moving in a time-dependent potential that is a super- 
position of analytic potentials of the individual galaxies. Such 
an approach maps the gravitational potential with high spatial 
resolution for low CPU costs due to the linear CPU scaling with 
the number of particles. However, the self-gravity of the stel- 
lar systems is not considered directly. That is, the orbital decay 
of the Magellanic Clouds due to dynamical friction cannot be 
treated self-consistently in restricted N-body simulations, but 
has to be considered by (semi-)analytical approximative for- 
mulas. 

The first papers on the physical features of the interacting 
system of the LMC, the SMC and the Galaxy used 3 D re- 
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Fig. 1. Contour map of the observed HI integrated relative 
column density in the Magellanic System. Data by Briins et 
al. ( 2005 1 is projected on the plane of sky. Galactic coordinates 
are used. 
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Fig. 2. LSR radial velocity of the Magellanic Stream as a func- 
tion of Magellanic Longitude (see Wannier & Wrixon 1972 1. 
The observational data by Briins et al. ( 2005 1 is plotted. Strong 
HI emission observed for i/lsr ~ Okms"' comes from the 
MW. The map also shows the emission from the SMC, the 
galaxies NGC 300, NGC 55, and NGC 7793 from the Sculptor 
Group, and the Local Group galaxy WLM. 



stricted N-body simulations to investigate the tidal origin of 
the extended Magellanic structures. Lin & Lynden-Bell (I1977> 
pointed out the problem of the large parameter space of the 
LMC-SMC-MW interaction. To reduce the parameter space, 
they neglected both the SMC influence on the System and dy- 
namical friction within the MW halo, and showed that such 



configuration allows for the existence of a LMC trailing tidal 
stream. The interaction between the Clouds was analyzed by 
Fujimoto & Sofue (WW), who assume the LMC and SMC to 
form a gravitationally bound pair for several Gyr, moving in 
a flattened MW halo. They identified some LMC and SMC 
orbital paths leading to the creation of a tidal tail. Following 
studies by Murai& Fujimoto (1980 1984 1, Lin & Lynden- 
Bell (il982L Gardiner et al. ( fT994t . and Lin et al. (2995) ex- 
tended and developed test-particle models of the LMC-SMC- 
MW interaction. The Magellanic Stream was reproduced as a 
remnant of the LMC-SMC encounter that was mostly placed 
at the time of -2 Gyr. The matter torn off" was spread along the 
paths of the Clouds. The simulations also indicate that the ma- 
jor fraction of the Magellanic Stream gas stems from the SMC. 
The observed radial velocity profile of the Stream was modeled 
remarkably well. However, the smooth H I column density dis- 
tribution did not agree with observations indicating apparently 
clumpy Magellanic Stream structure. Test-particle models 
place matter at the Leading Arm area naturally (see the study on 
the origin of tidal tails and arms by Toomre & Toomre ( I1972» . 
but coiTespondence with observational data cannot be consid- 
ered sufficient. Gardiner & Noguchi (1996) devised a scheme 
of the Magellanic System interaction implementing the full N- 
body approach. The SMC was modeled by a self-gravitating 
sphere moving in the LMC and MW analytic potentials. It 
was shown that the evolution of the Magellanic Stream and 
the Leading Arm is dominated by tides, which supports the 
applicability of test-particle codes for the modeling of ex- 
tended Magellanic structures. Recently, the study by Connors 
et al. (ii2005< investigated the evolution of the Magellanic 
Stream as a process of tidal stripping of gas from the SMC. 
Their high-resolution N-body model of the Magellanic System 
based on ideas of Gardiner & Noguchi (1996 1 is compared to 
the data from the HIPASS survey. Involving pure gravitational 
interaction allowed for remarkably good reproduction of the 
Magellanic Stream LSR radial velocity profile. They were able 
to improve previous models of the Leading Arm. Similarly 
to the previous tidal scenarios, difficulties remain concerning 
overestimations of the H I column density toward the far tip of 
the Magellanic Stream. Connors et al. ( 2005 1 approximate both 
the LMC and the MW by rigid potentials and also do not study 
the influence of the non-spherical halo of the MW. 

Meurer et al. (1985 1 involved continuous ram pressure 
stripping into their simulation of the Magellanic System. 
This approach was followed later by Sofue k\99Ai . who ne- 
glected the presence of the SMC, however. The Magellanic 
Stream was formed of the gas stripped from outer regions 
of the Clouds due to collisions with the MW extended ion- 



ized disk. Heller &Rohlfs ( fT994b argue for a LMC-SMC col- 
lision 500 Myr ago that established the MB. Later, gas dis- 
tributed to the inter-cloud region was stripped off' by ram pres- 
sure as the Clouds moved through a hot MW halo. Generally, 
continuous ram pressure stripping models succeeded in re- 
producing the decrease of the Magellanic Stream HI col- 
umn density towards the far tip of the Stream. However, they 
are unable to explain the evidence of gaseous clumps in the 
Magellanic Stream. Gas stripping from the Clouds caused by 
isolated collisions in the MW halo was studied by Mathewson 
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et al. ( I1984> . The resulting gaseous trailing tail consisted of 
fragments, but such a method did not allow for reproduc- 
tion of the column density decrease along the Magellanic 
Stream. Recently, Bekki&Chiba (2005 1 applied a complex 
gas-dynamical model including star-formation to investigate 
the dynamical and chemical evolution of the LMC. They in- 
clude self-gravity and gas dynamics by means of sticky parti- 
cles, but they are also not complete: they assume a live LMC 
system, but the SMC and MW were treated by static spheri- 
cal potentials. Thus, dynamical friction of the LMC in the MW 
halo is only considered by an analytical formula and a possible 
flattening of the MW halo is not involved. Their model cannot 
investigate the possible SMC origin of the Magellanic Stream 
either Mastropietro et al. (|2005 1 introduced their model of the 
Magellanic System including hydrodynamics (SPH) and a full 
N-body description of gravity. They studied the interaction 
between the LMC and the MW. Similarly to Lin&Lynden- 
BeU (fT977^ and Sofue ( 1994 1, the presence of the SMC was 
not taken into account. It was shown that the Stream, which 
sufficiently reproduces the observed H I column density distri- 
bution, might have been created without an LMC-SMC inter- 
action. However, the history of the Leading Arm was not inves- 
tigated. 

In general, hydrodynamical models allow for better repro- 
duction of the Magellanic Stream H I column density profile 
than tidal schemes. However, they constantly fail to reproduce 
the Magellanic Stream radial velocity measurements and espe- 
cially the high-negative velocity tip of the Magellanic Stream. 
Both families of models suffer from serious difficulties when 
modeling the Leading Arm. 

To model the evolution of the Magellanic System, the ini- 
tial conditions and all parameters of their interaction have to be 
determined. Such a parameter space becomes quite extended. 
In the Magellanic System we have to deal with the orbital 
parameters and the orientation of the two Clouds, their inter- 
nal properties (like the extension of the disk), and the prop- 
erties of the MW potential. In total we have about 20 pa- 
rameters (the exact number depends on the adopted sophis- 
tication of the model). Previous studies of the Magellanic 
Clouds, however, argued for very similar evolutionary scenar- 
ios of the system (e.g., Lin & Lynden-Bell 1982 Gardiner et 
al. ri994 Bekki&Chiba i2005i) . These calculations are based 
on additional assumptions concerning the orbits, or the inter- 
nal structure and orientation of the Clouds, the potential of 
the MW (mass distribution and shape), the treatment of dy- 
namical friction in the Galactic halo, or the treatment of self- 
gravity and gas dynamics in the Magellanic Stream. Some of 
them can be motivated by additional constraints. For example, 
Lin & Lynden-BeU J1982> and Irwin et al. J1990t argue that 
the Clouds should have been gravitationally bound over the last 
several Gyr. However, in general, the uniqueness of the adopted 
models is unclear, because a systematic analysis of the entire 
parameter space is still missing. 

We explore the LMC-SMC-MW interaction parame- 
ter space that is compatible with the observations of the 
Magellanic System available to date. Regarding the dimension 
and size of the parameter space, a large number of the numeri- 
cal model runs have to be performed to test possible parameter 



combinations, no matter what kind of search technique is se- 
lected. In such a case, despite their physical reliability, full N- 
body models are of little use, due to their computational costs. 

3. Method 

To cope with the extended parameter space, neither a complete 
catalog of models nor a large number of computationally ex- 
pensive self-consistent simulations can be performed, both due 
to numerical costs. However, a new numerical approach based 
on evolutionary optimization methods combined with fast (ap- 
proximative) N-body integrators turned out to be a promising 
tool for such a task. In case of encounters between two galaxies, 
Wahde J1998> and Theis J 19991 1 showed that a combination of 
a genetic algorithm with restricted N-body simulations is able 
to reproduce the parameters of the interaction. 

Here we apply the GA search strategy with a restricted N- 
body code for the Magellanic System. In the following sections 
we describe first our N-body calculations and then we briefly 
explain the applied genetic algorithm. 

3.1. N-body simulations 

Our simulations were performed by test-particle codes simi- 
lar to the ones akeady applied to the Magellanic System by 
Murai &Fujimoto (19801 and Gardiner et al. J1994> . But as an 
extension of these previous papers we allow for a flattened MW 
halo potential and a more exact formula for dynamical friction 
taking anisotropic velocity distributions into account. 

For the galactic potentials we used the following de- 
scriptions: both LMC and SMC are represented by Plummer 
spheres. The potential of the DM halo of the MW is 
modeled by a flattened axisymmetric logarithmic potential 
(Binney & Tremaine 1987) 

OL = ^i^o'ln(^c'+/?' + ^). (1) 

In agreement with Helmi ( I2004> we set = 12kpc and vq - 
V2 ■ 13 1.5 km s"', and q describes the flattening of the MW 
halo potential. The corresponding flattening qp of the density 
distribution associated with the halo potential follows 

qj^ = ^ "^^^^ (R <K R^), (2) 
2 + 3/^2 V cy, w 

^p2^ ^2^2-^1 {R:i>R,). (3) 

The choice for the logarithmic potential was motivated by sev- 
eral reasons. First, it allows for the investigation of flattened 
configurations of the Galactic halo at low computational costs. 
Also the relatively small number of input parameters of Eq. ^ 
is an advantage concerning the performance and speed of the 
numerical model. Finally, the logarithmic potential was em- 
ployed in the recent study by Helmi ( 200^ that used the similar 
method of dwarf galaxy streams to investigate the MW DM 
halo, and so the application of the same formula allows for 
comparison of our results. 
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Dynamical friction causes the orbital decay of the 
Magellanic Clouds. We adopted the analytic dynamical fric- 
tion formula by Binney J1977I I. In contrast to the commonly 
used expression by Chandrasekhar (1943), it allows for an 
anisotropic velocity distribution. By comparison with N-body 
simulations of sinking satellites, Penarrubia et al. (2004i 
showed that Binney's solution is a significant improvement 
over the standard approach with Chandrasekhar's formula. 

Finally, we get the following equations of motion of the 
Clouds: 



dv] 



LMC 



df 



-(VcDl + VOsMc) + 



'DF 

rriLMC ' 



df 



msMC 



(4) 



(5) 



where Olmc, are the LMC, SMC Plummer potentials, 

and Fdf is the dynamical friction force exerted on the Clouds 
as they move in the MW DM halo. 

Our simulations were performed with the total number of 
10000 test-particles equally distributed to both Clouds. We 
start the simulation with test-particles in a disk-like configura- 
tion with an exponential particle density profile, and compute 
the evolution of the test-particle distribution up to the present 
time. Initial conditions for the starting point of the evolution 
of the System at the time -4 Gyr were obtained by the standard 
backward integration of equations of motion (see, e.g., Murai & 
Fujimoto 1980 Gardiner et al. 1994). Basically, the choice for 
the starting epoch of this study originates in the fact that the 
MW, LMC, and SMC form an isolated system in our model. 
Such an assumption was very common in previous papers on 
the Magellanic System (e.g., Murai & Fujimoto 1980, Gardiner 



et al. 119941 Gardiner & Noguchi 1996) and was motivated by 
insufficient kinematic information about the Local Group, that 
did not allow for the estimation of the influence of its mem- 
bers other than MW on the evolution of the Magellanic System. 
Our detailed analysis of the orbits of LMC and SMC showed 
that the galactocentric distance of either of the Clouds did not 
exceed 300 kpc within the last 4 Gyr Investigation of the kine- 
matic history of the Local Group by Byrd et al. ( 1994i indicates 
that the restriction of the maximal galactocentric radius of the 
Magellanic System to /^max ~ 300 kpc when T > -4 Gyr lets us 
consider the orbital motion of the Clouds to be gravitationally 
dominated by the MW. 

3.2. Genetic algorithm search 

Genetic algorithms can be used to solve optimization problems 
like a search in an extended parameter space. The basic con- 
cept of GA optimization is to interpret the natural evolution 
of a population of individuals as an optimization process, i.e., 
an increasing adaptation of a population to given conditions. 
In our case the conditions are to match numerical models to 
the observations. Each single point in parameter space uniquely 
defines one interaction scenario that can be compared with the 
observations (after the N-body simulation is performed). The 
quality of each point in parameter space (or the corresponding 
N-body model) can be characterized by the value of a fitness 



function (FF), which is constructed to become larger, the bet- 
ter the model matches the observations. A population consists 
of a set of points in parameter space (individuals). Each single 
parameter of an individual corresponds to a gene. A genetic 
algorithm recombines the genes of the individuals in differ- 
ent reproduction steps: First two individuals (parents) are se- 
lected with a probability growing with their fitness. Then the 
genes of the parents are recombined by application of repro- 
duction operators mimicking cross-over and mutation. Often 
the reproduction is done for all members of a population. Then 
the newly created population corresponds to a next generation. 
The reproduction steps are then repeated until a given num- 
ber of generations is calculated or a sufficient convergence is 
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Fig. 3. Contour map of the observed H 1 integrated column den- 
sity. The original data cube by Briins et al. (2005) is modified 
by frequency filtering (see Appendix^ and integrated over all 
radial velocity channels. The contour map is projected on the 
plane of sky. Galactic coordinates are used. 

reached. More details about genetic algorithms can be found 



in Holland (il975i> or Goldberg ( I1989> . An application to inter- 
acting galaxies is described in Wahde |1998 1, Theis ( 1999 1, or 
Theis &Kohle (20011. 

Obviously, the evolutionary search for the optimal solution 
can be treated as a process of maximizing the fitness of individ- 
uals. The GA looks for maxima of the function assigning indi- 
viduals their fitness values. It should be noted that the problem 
we want to solve only enters via the FF. Therefore, the choice 
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of the FF is essential for the performance and the answer given 
by a GA. 

For our calculations we used a FF consisting of three 
different parts corresponding to three different comparisons. 
These three FFs measure the quality of the numerical mod- 
els for different aspects of the given HI data cube. The origi- 
nal 3D H I data cube together with its version modified for the 
purpose of an efficient GA search are visualized in Fig. IB. II 
Two of the comparisons deal with the whole data cube: FF2 
denotes the rough occupation of cells in the data cube and is 
a measure for the agreement of the structural shape in the data 
cube. FFt, compares the individual intensities in each cell of 
the data cube. The FFt, definition basically follows the FF in- 
troduced in Theis <1999> or Theis & Kohle J2001I I. who found 
it to be an efficient GA driver for the galactic interaction prob- 
lems they studied. However, if the fraction of the total volume 
of the system's data cube that is occupied by the structures of 
special interest is small (< 10% in the case of the Magellanic 
Stream and the Leading Arm), implementation of a structural 
search in the data (FF2) preceding the fine comparison between 
modeled and observed data significantly improves the perfor- 
mance of the GA. Finally, FFi is introduced to take into ac- 
count the velocity profile, i.e., the important constraint of the 
minimum velocity in the stream. All these three quality mea- 
sures are combined to yield the final fitness of a model. Details 
can be found in AppendixiBl 



Table 1. Parameters of the MW-LMC-SMC interaction. 



Param. 



Value 



J-LMcLkpc] 



'■sMc[kpc] 



<-1.5,-0.5> ^ 
(-41.0, -40.0) 
(-27.1, -26.1) , 

<13.1, 14.1) ^ 
(-34.8,-33.8) 
(-40.3, -39.3) 



Current galactocentric 
position vectors 



WLMc[kms 



WsMc[kms 



<-3,85> 
<-23I,-169> 
(132, 206) 

(-112,232) 
(-346, -2) 
(45,301) 



Current velocity 



vectors 



mLMc[10^'Mo] 


(15.0, 25.0) 


Masses 


msMcLlO'Mo] 


(1.5,2.5) 




fLMc[kpc] 


(2.5,3.5) 


Plummer sphere 


6sMc[kpc] 
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3.3. The GA parameter space 

In this paragraph we introduce the parameters of the 
Magellanic System interaction that were subject to our GA 
search. The parameters involve the initial conditions of the 
LMC and SMC motion, their total masses, parameters of mass 
distribution, particle disc radii, and orientation angles, and also 
the MW dark matter halo potential flattening parameter. Tabled 
summarizes the parameters of the interaction and introduces 
their current values. Models are described in a right-handed 
Cartesian coordinate system with the origin in the Galactic cen- 
ter. This system is considered to be an inertial frame because 
we assume that njMw » 'Mlmc, ^msmc, where niMw is the mass 
of MW, and wlmc and otsmc are the masses of LMC and SMC, 
respectively. Therefore, the center of mass of the system may 
be placed at the Galactic center. We assume the present posi- 
tion vector of the Sun Rq - (-8.5, 0.0, 0.0) kpc. In the follow- 
ing paragraphs we will discuss the parameters of the LMC- 
SMC-MW interaction and the determination of their ranges. 

The nature and distribution of dark matter in the Galaxy has 
been subject to intense research and a large number of mod- 
els have been proposed. We probe the DM matter distribution 
by investigating the redistribution of matter in the Magellanic 
System due to the MW-LMC-SMC interaction, paying special 
attention to the features of the Magellanic Stream. That is sim- 
ilar to the method applied by Helmi ( 2004j . who studied kine- 
matic properties of the Sagittarius stream. To enable compari- 
son with the results by Helmi (2004), we also adopted the axi- 
ally symmetric logarithmic halo model of MW (Eq.^ and the 
same values of the halo structural parameters Rc, and vq with a 



similar range of studied values of the flattening q (see TabIe[Q. 
We extended the range of q values tested by Helmi ( I2004> to the 
lower limit of ^ = 0.74, which is the minimal value acceptable 
for the model of a logarithmic halo (for a detailed explanation 
see Binney & Tremaine ( 1987)). For every value of q, a time- 
consuming calculation of parameters of dynamical friction is 
required (see Appendix|A|i. To reduce the computational diffi- 
culties, the flattening q was treated as a discrete value with a 
step of Aq = 0.02, and the parameters of dynamical friction 
were tabulated. The upper limit of ^ = 1.20 was selected to 
enable testing of prolate halo configurations. Extension of the 
halo flattening range to higher values was not considered due 
to the computational performance limits of our numerical code. 

The estimated ranges of the values of the remaining pa- 
rameters are based on our observational knowledge in the 
Magellanic System. Galactocentric position vectors Tlmc and 
rsMC agree with the LMC and SMC distance moduli measure- 
ments given by Van den Bergh ( 2000 1, who derived the mean 
valu es of dist ance moduH (m - M)q = 18.50 + 0.05 for LMC 
and (m - M)q - 18.85 ± 0.1 for SMC, corresponding to the 
heliocentric distances of (50.1 + 1.2) kpc and (58.9 ± 2.6) kpc, 
respectively. Only 2 of the 6 components of the LMC and SMC 
position vectors enter the GA search as free parameters be- 
cause the rest of them are determined by the projected posi- 
tion of the Clouds on the plane of sky, that is Iimc - 280° 27', 
W = -32° 53' and Zsmc = 302° 47', W = -44° 18'. 

Previous studies by Murai & Fujimoto J 19801 1 and Gardiner 
et al. Cl994i) found that the correct choice of the spatial 
velocities is crucial for reproducing the observed HI struc- 
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tures. Kroupa & Bastian < 19971 1 derived proper motions of 
LMC and SMC stars from an analysis of HIPPARCOS such 
that t/LMC = (+41+44,-200±31,+169±37)kms-' and 
vsMC = (+60+ 172,-174+ 172, +173 ± 128) km s"'. The 
large uncertainties in actual values of the LMC and SMC ve- 
locity vectors originate in the fact that the measured transverse 
velocities suffer from large rms errors, which is connected 
to the large distance of LMC and SMC in combination with 
the limited precision of proper motions in the HIPPARCOS 
catalog. To derive the actual initial conditions at the starting 
time of the simulation from the current positions and velocities 
of the Clouds, we employed the backward integration of 
equations of motion. 

Current total masses otlmc and msMC follow estimates by 
Van den Bergh (2000). In general, masses of the Clouds are 
functions of time and evolve due to the LMC-SMC exchange 
of matter, and as a consequence of the interaction between 
the Clouds and MW. Our test-particle model does not allow 
for a reasonable treatment of a time-dependent mass-loss. 
Therefore, masses of the Clouds are considered constant in 
time, and their initial values at the starting epoch of our simula- 
tions are approximated by the current LMC and SMC masses. 
Regarding the large range of the LMC and SMC mass esti- 
mates available (for details see Van den Bergh (2000)), /wlmc 
and msMC are also treated as free input parameters of our model 
that become subjects to the GA optimization. Their ranges that 
we investigated can be found in Tabled 

Scale radii of the LMC and SMC Plummer potentials e^MC 
and esMC are input parameters of the model describing the ra- 
dial mass distribution in the Clouds. The study by Gardiner et 
al. Il994i used the values fLMC = 3kpc and esMC - 2kpc. 
To investigate the influence of this parameter on the evolution 
of the Magellanic System, the Plummer radii were involved in 
the GA search and their values were varied within the ranges 
of the width of 1 kpc, including the estimates by Gardiner et 
al. ( fl994.) (see Tab leffli. 

Gardiner et al. ( I1994> analyzed the H I surface contour map 
of the Clouds to estimate the initial LMC and SMC disk radii 
entering their model of the Magellanic System as initial con- 
ditions. Regarding the absence of a clearly defined disk of the 
SMC, and possible significant mass redistribution in the Clouds 
during their evolution, the results require careful treatment and 
further verification. We varied the current estimates of disk 



radii '"lmc '"smc ^^^^ parameters within the ranges in- 
troduced in Tabled containing the values derived by Gardiner 
et al. ( I1994> . and used them as initial values at the starting point 
of our calculations (T - -4Gyr). 

The orientation of the disks is described by two angles © 



and O defined by Gardiner & Noguchi ( I1996t . Several observa- 
tional determinations of the LMC disk plane orientation were 
collected by Lin et al. (,1995J . Its sense of rotation is assumed 
to be clockwise (Lin et al. 1 19951 Kroupa & Bastian 1997). The 
position angle of the bar structure in the SMC was used by 
Gardiner & Noguchi ( 1996) to investigate the current disk ori- 
entation. Their results allow for wide ranges of the LMC and 
SMC disk orientation angles (see Tabled and so we inves- 
tigated © and C) by the GA search method, too. Similarly to 
Gardiner & Noguchi fl996j . we use the current LMC and SMC 



disk orientation angles (Table [TJ to approximate their initial 
values at r = -4 Gyr. 

4. Results 

We try to reproduce as closely as possible the column density 
and velocity distribution of H I in the Magellanic Stream and 
in the Leading Arm. The influence of actual distances to the 
LMC and SMC and of their present space velocity vectors is 
considered together with their masses and the past sizes and 
space orientations of the orginal disks. Here, we give the results 
of the search in the parameter space with the GA using the 3- 
component FF defined by Eg. IB. II In principle, the GA is able 
to find the global extreme of the FF if enough time is allowed 
for the evolution of the explored system (see Holland ri975> . 
However, it may be very time-consuming to identify the single 
best fit due to a slow convergence of the FF. Therefore, to keep 
the computational cost reasonable, the maximum number of 
120 GA generations to go through was defined. 

To explore the FF of our system, we collected 123 GA 
fits of the Magellanic System resulting from repeated runs of 
our GA optimizer. Typically, identification of a single GA fit 
requires x 10'^ runs of the numerical model. Thus, due to the 
application of GA we were able to search the extended pa- 
rameter space of the interaction and discover the most suc- 
cessful models of the System over the entire parameter space 
by testing ^ lO*" parameter combinations. In the case of our 
20-dimensional parameter space, simple exploration of every 
possible combination of parameters even on a sparse grid of, 
e.g., 10 nodes per dimension, means 10^*^ runs of the model. 
Such a comparison clearly shows the necessity of using opti- 
mization techniques and demonstrates computational efficiency 
ofGA. 




0.8 
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1.0 
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1.2 



Fig. 4. Maximum values of fitness as a function of the MW dark 
matter halo flattening q. The plot depicts the fitness of the best 
GA fit of the Magellanic System that was found for each of 
the MW dark matter halo flattening values that entered the GA 
search. The function is also approximated by its least-square 
polynomial fit. The values of q delimiting the model groups A, 
B, and C (see Tabled are emphasized by dotted Unes. 

The fitness distribution for diff'erent values of the halo flat- 
tening parameter q is shown in Fig. |4] For every value of q, 
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the model of the highest fitness is selected and its fitness value 
is plotted. Figure 0] indicates that better agreement between 
the models and observational data was generally achieved 
for oblate halo configurations rather than for nearly spheri- 
cal or prolate halos. The relative difference between the worst 
(F;fMiN = 0.450;^ = 1.14) and the best (FituAX - 0.496; = 
0.84) model shown in Fig. His AFit - 0.09. It reflects the 
fact that each of the GA fits contains a trailing tail (Magellanic 
Stream) and a leading stream (Leading Arm), but that there are 
fine differences between the resulting distributions of matter. 
One may note that the GA optimizer did not discover a model 
of fitness over 0.5 (the maximum reachable fitness value is 1.0 
- see Appendix IbI. It is caused either by insufficient volume 
of the studied parameter space of the interaction, or by simpli- 
fication of physical processes in our model (see Sec.|6j, or by 
a combination of both reasons. That establishes an interesting 
problem and should become a subject to future studies. 

We want to discuss our results with respect to the shape 
of the MW halo. Thus, all the GA fits are divided into three 



Table 2. Three major groups according to the halo flattening q 



Group 


A 


B 


C 




0.74 <q< 0.92 


0.94 < ^ < 1.06 


1.08 <q< 1.20 


N, 


101 


10 


12 



groups according to the halo flattening (see Table ^ to show 
differences or common features of models for oblate, nearly 
spherical, and prolate halo configurations. The actual borders 
between the groups A, B, and C were selected by definition 
and, so that the number of models in each of the groups allows 
for statistical treatment of the LMC and SMC orbital features 
and particle redistribution that will be introduced in Sect. 14.11 

4.1. Evolution 

Close encounters of interacting galaxies often lead to substan- 
tial disruption of their particle disks, forming subsequently 
tidal arms and tails (Toomre & Toomre ll972> . Regarding that, 
the time dependence of the relative distance of the interacting 
pair is an interesting source of information about the system. 

First of all, we examine the time distribution of the min- 
ima of the LMC-SMC relative distance. For each of the model 
groups mentioned in Tabled we calculate the relative number 
of GA solutions having a minimum of the LMC-SMC relative 
distance within a given interval of 500 Myr Figure |5] shows 
such a distribution of fits for the total time interval of our sim- 
ulations, which is 4.0 Gyr. The local maxima of the time dis- 
tribution of the LMC-SMC distance minima are within the in- 
tervals (-3.0, -2.0) Gyr and (-0.5, 0.0) Gyr. For prolate halos 
(q > 1.08) there is no LMC-SMC distance minimum between 
-2.0 Gyr and -0.5 Gyr. 

Subsequently, the mean values of the LMC-SMC distance 
minima are calculated for each of the time intervals defined 
above. Comparison of the results for oblate, nearly spherical, 
and prolate DM halo configuration is available in Fig. |6l It 
was found that close (Ar « lOkpc) LMC-SMC encounters 



do not occur in models with either oblate or prolate halos. If 
the MW DM halo shape is nearly spherical, disruption of the 
LMC and SMC initial particle distribution leading to creation 
of the observed H I structures typically occurred due to strong 
LMC-SMC interaction. 

Another point of interest is the time dependence of the 
LMC and SMC test-particle redistribution during the evolu- 
tionary process. Figure offers the relative number of test- 
particles strongly disturbed, i.e., particles that reached the min- 
imal distance of twice the original radii of their circular or- 
bits around the LMC and SMC centers, respectively, and by 
the LMC-SMC-MW interaction in the defined time-intervals. 
Comparison between the plots in Figs. |6l and shows that 
encounters of the Clouds are followed by delayed raise of 
the number of particles shifted to different orbits, typically. 
Another such events are induced by the interaction of the 
Clouds and MW. Disruption of the LMC and SMC disks trig- 
gers formation of the extended structures of the Magellanic 
System. Particles are assigned new orbits in the superimposed 
gravitational potential of the LMC, SMC, and MW, and spread 
along the orbital paths of the Clouds. Our study shows that the 
formation of the Magellanic Stream and other observed H I fea- 
tures did not begin earlier than 2.5 Gyr ago for model groups 
A and B (see Fig.0. Prolate halos (group C) allow for a mass 
redistribution in the system that started at T < -3.5 Gyr 

4.2. Representative models 

Here, we describe the models of highest fitness selected from 
each of the groups A, B, and C. All of them are typical repre- 
sentatives of their model groups and we discuss their features 



Table 3. Parameters of the best models in separate q categories. 
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Fig, 5. Relative numbers of the IVIagellanic System GA models with LMC-SIVIC relative distance minima at a given time interval 
for model groups A (left), B (middle), and C. The counts were made for 8 time intervals of 500 Myr covering the entire Magellanic 
evolution period of 4 Gyr investigated in our study. 
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Fig. 6. Mean values of the LMC-SMC distance minima for model groups A (left), B (middle), and C (Table |2j and for 8 time 
intervals of 500 Myr. 
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Fig. 7. Relative number of the LMC and SMC test-particles strongly disturbed due to the Magellanic System interaction. Counts 
within 8 time intervals of 500 Myr and for model groups A (left), B (middle), and C are plotted. 



with respect to the HI observational data. Table |3l summarizes 
the parameter values of the models. 

4.2.1. Group A 

The best model of the Magellanic System with an oblate MW 
halo is introduced in this section (model A). Figure |S] depicts 
the time variation of the LMC and SMC galactocentric dis- 
tances together with the LMC-SMC separation for the last 
4 Gyr. The Clouds move on very different orbits. The apogalac- 
tic distance of LMC decreases systematically during the evolu- 
tionary period, which clearly reflects the effect of dynamical 
friction. There is a gap between the periods of subsequent peri- 
galactic approaches of the Clouds. While the last two perigalac- 
tica of LMC are separated by » 2.3 Gyr, it is not over » L5 Gyr 
in the case of the SMC. Filled parts of the plot in Fig.|8]indicate 



that the Magellanic Clouds have reached the state of a gravita- 
tionally bound system during the last 4 Gyr. We define grav- 
itational binding by the sum of the relative kinetic and grav- 
itational potential energy of LMC and SMC. The Clouds are 
bound when the total energy is negative. Specifically, LMC and 
SMC have been forming a bound couple since T - -L06 Gyr. 
Nevertheless, the total lifetime of a bound LMC-SMC pair did 
not exceed 40 % of the entire evolutionary period we studied. 

Comparison between Figs. |H1 and |5] allows for conclusions 
about major events that initialized the redistribution of the 
LMC and SMC particles. The most significant change of the 
initial distribution of particles occurred as a result of the LMC- 
MW approach at T a; -2.4 Gyr, and the preceding LMC-SMC 
encounter T ^ -2.5 Gyr. Later on, the particle redistribution 
continued due to tidal stripping by the MW. 
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FigureQ^shows the modeled distribution of integrated HI 
column density in the System. To enable comparison with the 
observed H I distribution, we plotted a normalized H 1 column 
density map. The technique used to convert a test-particle dis- 
tribution into a smooth map of column densities is described 
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Fig. 8. Orbital evolution of the Clouds for the best GA fit from 
the model group A. The plot corresponds to an oblate halo of 
flattening q - 0.84. Time dependence of the LMC (dashed line) 
and SMC (dotted line) galactocentric distances, and the LMC- 
SMC relative distance are plotted above. Plot areas with grey 
filling mark the time intervals when the Clouds were gravita- 
tionally bound to each other. 

in AppendixiBl Mass distribution of H I extends beyond the far 
tip of the observed Magellanic Stream (Fig.|3} in the model A. 
H I column density peaks can be found in Fig. |3l at the posi- 
tions / = 300°, b = -65° and / = 45°, b = -82°. The model 
A places local density maxima of H I close to those observed 
ones (i.e., relative angular distance is ^ 10°) to approximate po- 
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Fig. 9. Model A. Relative number of LMC/SMC test-particles 
strongly disturbed due to the interaction in the Magellanic 
System. Counts within 8 time intervals of 500 Myr are plotted. 



sitions I = 325°, b = -70° and / = 70°, b = -70°, respectively. 
Note also the low-density distribution of matter spread along 
the great circle of / = 270° (Fig.llO>. The matter emanates from 
LMC near the position of the Interface Region identified by 
Briins et al. l2005l ("Fig.mi. In general, the model overestimates 
the amount of matter in the Magellanic Stream. The Leading 
Arm consists only of LMC particles in this scenario. The mod- 
eled matter distribution covers a larger area of the plane of 
sky than what is observed. However, this is a common prob- 



315° 



270° 



+45° 




1.000 


■ 


+ oU 










1 1 


0.875 




+ 15° 












0° 


, +- ' 1 " ^ 


0.750 




-15° 




0.625 




-30° 










0.500 




-45° 




0.375 














' \ ,^ ~^ r 


0.250 




-45° 












0.125 




-30° 












0.000 





90° 



Fig. 10. Model A. Contour map of the modeled H I integrated 
relative column density. Data is projected on the sky plane. 
Galactic coordinates are used. 

lem of previous test-particle models of the Magellanic System 
(see, e.g., Murai & Fujimoto 1 980. Gardiner et al. 1994 1 and is 
likely caused by simplifications in the treatment of the physi- 
cal processes. But also in general, successful reproduction of 
the Leading Arm has been a difficult task for all the models 
introduced so far 

The LSR radial velocity profile of the Magellanic Stream 
for the model A is shown in Fig. ^2 The model reproduces 
the LSR radial velocity of the Stream matter as an almost lin- 
ear function of Magellanic Longitude with the high-negative 
velocity tip reaching -400 km s"'. Such features are in agree- 
ment with observations (see Fig. |2j. In contrast to Gardiner 
et al. (19941, the Magellanic Stream consists of both LMC 
and SMC particles. Figure[TT]denotes that the LMC and SMC 
Stream components cover different ranges of LSR radial ve- 
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Fig. 11. Model A. LSR radial velocity profile of the Magellanic 
Stream. 



locities. The Stream component of the SMC origin does not 
extend to LSR radial velocities below the limit of -200 km s ' . 
The major fraction of the LMC particles resides in the LSR 
radial velocity range from -400 km s"' to -200 km s"'. 



4.2.2. Group B 

The best model of the group B (model B) corresponds to the 
MW DM halo flattening value q - 1.02. The initial condition 
set for the model B is listed in Table|3] The galactocentric dis- 
tance of the Clouds and their spatial separation as functions 
of time are plotted in Fig.[^ Continuous decrease of the LMC 
and SMC galactocentric distances due to the dynamical friction 
is apparent for both LMC and SMC. A very close encounter of 
the Clouds with the relative distance Ar a; lOkpc occurred at 
r a -2.2 Gyr At similar moments ofT a -2.3 Gyr (LMC) and 
T a -2.1 Gyr (SMC), the Clouds also reached perigalactica of 
their orbits. In general, both Clouds have been moving on orbits 
showing similar time dependence of their galactocentric dis- 
tances, as indicated by Fig. El Nevertheless, the position vec- 
tors of the Clouds evolved in significantly different ways. As a 
consequence, the spatial separation of the Clouds varied within 
a wide range of values from Ar ^ 10 kpc to Ar ^ 250 kpc. The 
Clouds have formed a gravitationally bound couple three times 
within the last 4 Gyr, and the total duration of such periods was 
1 .7 Gyr. Currently, LMC and SMC are gravitationally bound in 
the model B. 

The LMC-SMC encounter at T » -2.2 Gyr caused a dis- 
tortion of the original particle disks of the Clouds. More than 
25% of the total number of the LMC and SMC particles 
were moved to significantly different orbits (for a definition 
see Sect. 14. Il l within the interval of 1 Gyr after the encounter 
(see Fig.ll3>. The following evolution of the particle distribu- 
tion formed extended structures depicted in Fig.^^ There are 
two spatially separated components present in the modeled tail. 
The H 1 column density distribution map for the model B (see 
Fig. I14l i shows a densely populated trailing stream parallel to 
the great circle of I - 1357315°. It consists of the SMC par- 



500 



400 



o 



^ 300 



O 

c 

CO 

« 

b 



B 200 



100 







q = 1 .02 

LMC-MW 
SMC-MW 
LMC-SMC 




-4 -3 -2 

Time [Gyr] 







Fig. 12. Orbital evolution of the Clouds for the best GA fit 
from the model group B. The plot corresponds to a nearly 
spherical halo of flattening q - 1.02. Time dependence of the 
LMC (dashed line) and SMC (dotted line) galactocentric dis- 
tances, and the LMC-SMC relative distance are plotted above. 
Grey filled areas show time intervals when the LMC and SMC 
formed a gravitationally bound couple. 
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Fig. 13. Model B. Relative number of LMC/SMC test-particles 
strongly disturbed due to the interaction in the Magellanic 
System. Counts within 8 time intervals of 500 Myr are plotted. 



tides torn off from the initial disk ^ 2 Gyr ago. Its far end is 
projected to the plane of sky close to the tip of the Magellanic 
Stream. However, the modeled increase of the column density 
of matter toward the far end of the tail is a substantial draw- 
back of scenario B. The stream extends into the SMC leading 
arm located between / = 290°, ^ = -15° and / = 290°, b = 45°. 
The second component of the particle tail is of LMC origin 
and is spread over the position of the observed low-density gas 
distribution centered at / = 80°, b = -50° (Figs.[T]or|3}. 

The most significant structure at the leading side of the 
Magellanic System is the SMC stream introduced in the pre- 
vious paragraph. Comparison between Figs.01andll4lindicates 
that neither the projected position nor the integrated H I den- 
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sity distribution of the stream is in agreement with the observed 
Leading Arm. There was also a structure emanating from the 
leading edge of the LMC identified at approximate position 
/ = 210°, b = -15° (see Fig. [HJ. Regarding the HI data by 



LSR radial velocity profile of the trailing tail of the model B 
does not extend over the limit of^LSR ~ -350 km. s"' (Fig. ll5> . 
However, following the HI data, the far tip of the Stream 
should reach the LSR radial velocity ulsr ~ -400 km. s ' at the 
Magellanic Longitude a; -100°. 

4.2.3. Group C 

Our last model group C assumes the presence of prolate MW 
DM halos. The best GA fit of the System (model C) is intro- 
duced in Table|3lreviewing its initial condition set. Concerning 
orbital motion of the Clouds, there is significant diff'erence be- 
tween the LMC and SMC periods of perigalactic approaches 
during the last 4Gyr While the last period of the LMC ex- 
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Fig. 16. Orbital evolution of the Clouds for the best GA fit from 
the model group C. The plot corresponds to a prolate halo of 
flattening q = 1.16. Time dependence of the LMC (dashed 
line) and SMC (dotted line) galactocentric distances, and the 
LMC-SMC relative distance are plotted above. Periods when 
the Clouds formed a gravitationally bound couple are marked 
by grey filling. 

ceeds 2.5 Gyr, the SMC orbital cycle is shorter than 1.5 Gyr. 
The relative distance of the Clouds remains over 70kpc for 
T < -0.4 Gyr. They became a gravitationally bound couple 
0.6 Gyr ago and this binding has not been disrupted (Fig.ll6>. 

Changes to the original LMC and SMC particle disks oc- 
curred especially due to the LMC-MW and SMC-MW en- 
counters at r < -2.0 Gyr. Comparison between Figs.[T6land[T7l 
demonstrates the significance of different encounter events for 
particle redistribution. Note that the rise of the number of dis- 
turbed particles is delayed with respect to the corresponding 
disturbing event. Subsequently, the evolution of particle struc- 
tures continued under the influence of tidal stripping by the 
gravitational field of MW. The current distribution of matter in 
the model C is plotted in the form of a 2-D map in Fig.^] The 
projection of the modeled trailing stream indicates that it occu- 
pies larger area of the data cube than the observed Magellanic 
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Fig. 17. Model C. Relative number of LMC/SMC test-particles 
strongly disturbed due to the interaction in the Magellanic 
System. Counts within 8 time intervals of 500 Myr are plotted. 
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Fig. 18. Model C. Contour map of the modeled H I integrated 
relative column density. Data is projected on the sky plane. 
Galactic coordinates are used. 



Stream (compare Figs. |3 and I18> . According to Fig. |3 the 
Magellanic Stream shows HI density peaks at I - 300°, 
b = -65° and / = 45°, b = -82°. Our model C expects 
two local maxima of H I integrated column density in the tail. 
Their positions are shifted by ^ 20° relatively to the peaks in 
Fig. 13 Additional comparison between the model and obser- 
vations reveals that the model C overestimates the integrated 



column densities of H I in the Magellanic Stream. The matter 
located at the leading side of the Magellanic System is of SMC 
origin only. Similarly to the case of the trailing tail, the mod- 
eled amount of matter exceeds observational estimates for the 
Leading Arm. 

Figure offers the LSR radial velocity profile of the 
Magellanic Stream in our model C. The measured minimum 
of the LSR radial velocity is x -400 km s ' . The high negative 
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LSR radial velocity at the far tip of the modeled Magellanic 
Stream does not exceed » -300 km s ', however The observed 
H I emission intensity decreases towards the high negative ve- 
locity tip, which was not well reproduced by the model C. 

5. Summary of the GA models 

5.1. Orbits of the Magellanic Clouds 

Exploration of the orbital motion of the Clouds shows the sim- 
ilarity of the GA fits for oblate and prolate halos (models A 
and C). No close (r < lOkpc) LMC-SMC encounters occurred 
for either of the models A and C. It is also notable that in the 
models with aspherical halos, the SMC period of perigalactic 
approaches is significantly shorter than the period of the LMC 
and that the SMC remains closer than lOOkpc to the MW cen- 
ter during the last 4 Gyr. When the MW DM flattening q ^ 1.0, 
the LMC and SMC orbital cycle lengths were comparable for 
the model B. Independently of the MW halo shape, the LMC 
and SMC are currently forming a gravitationally bound cou- 
ple in our models. However, the Clouds cannot be considered 
bound to each other during the entire period of the last 4 Gyr. 
This is in contrast with Murai & Fujimoto ( 1980 1, Gardiner et 
al. (lT994t . or Gardiner & Noguchi ( 1996 1, who argue that the 
LMC and SMC moving in the spherical halo have formed a 
gravitationally bound pair for at least several Gyr to allow for 
sufficient matter redistribution. We show that the structural res- 
olution adopted by the above-cited studies to make a compar- 
ison between models and observations does not allow for such 
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constraints of the orbital history of the Clouds. The GA search 
employed a 3-level detailed evaluation of the modeled HI 
distribution with respect to high-resolution observational data 
(Briins et al. 2005 1, which introduced significant improvement 
of previous approaches to compare observations and models. 
Nevertheless, there is still no clear indication that continuous 
gravitational binding of the Clouds covering the entire evolu- 
tionary period is necessary for successful reproduction of the 
observed data. 

5.2. Origin of the matter in the Stream 

In our best model both SMC and LMC particles were present 
in the trailing stream. This is a common feature of the scenarios 
that were investigated. In general, the fraction of H I gas orig- 
inating at the SMC exceeds the fraction of LMC matter in the 
Stream. 

Following the models A and B, the formation of the 
Magellanic Stream did not start earlier than 2.5 Gyr ago. In 
the case of the models C, the age of the Stream is » 3.5 Gyr. 
The estimates by the models A and B are close to the epoch 
when a massive star formation burst in LMC began (Van den 
Berghl^nnn^. It is also true for the model C. However, the dis- 
turbing event that triggered the formation of the Magellanic 
Stream in the model C occurred at a time that does not substan- 
tially differ from the epoch T = -4 Gyr when our simulation 
started. In such a case we have to assume that the tidal interac- 
tion responsible for the creation of the Stream was very likely 
invalidating the natural condition of sufficient isolation of the 
Clouds at the starting point of the simulation. Such a condi- 
tion is necessary to set up the initial particle disks around LMC 
and SMC when the particles reside on circular orbits. Thus, 
the statements on the actual age of the Stream in the model C 
should be made with care. Regardless, the large uncertainty of 
the estimated beginning of the LMC star burst (see Van den 
Bergh 2000 1 requires that consideration is given to the possi- 
bility that for case A, B, or C the starting epochs of the Stream 
creation and of the start formation burst in the LMC may over- 
lap. 

The previous paragraph indicates that it is very likely that 
matter forming the Stream comes from the Magellanic Clouds 
containing stars, and we necessarily face the observational fact 
that there is no stellar content in the Magellanic Stream. The 
models for aspherical halos indicate that the matter coming to 
the Magellanic Stream from the LMC originates in the outer 
regions of its initial particle disk, while no matter was torn off 
from the inner disk of radius rjisk ~ 5 kpc that was the dominant 
region of star formation in the LMC. It is due to the absence of 
close encounters in the Magellanic System. 

In contrast to the models A and C, a dramatic encounter 
event between the Clouds occuiTed in model B at s; -2.2 Gyr, 
when the internal structure of both disks was altered and the 
matter from central areas of the LMC disk was also transported 
to the Magellanic Stream. In such a case we expect a certain 
fraction of the matter of the Magellanic Stream to be in the 
form of stars, which is, however, not supported by observations. 



5.3. Structure of the Stream 

Mathewson et al. <1977> observationally mapped the 
Magellanic Stream and discovered its clumpy structure 
consisting of six major HI clouds named MSI- VI. Recently, 
a more sensitive high-resolution H I survey of the Magellanic 
System by Briins et al. J2005t showed that the above men- 
tioned fragments of the Stream have to be considered density 
peaks of an otherwise smooth distribution of neutral hydrogen 
with a column density decreasing towards the high-negative 
radial velocity tip of the Magellanic Stream. Our models 
corresponding to aspherical MW halos (A, C) placed local 
density maxima of H I close to the South Galactic pole. That 
result is supported by observations by Briins et al. (?2005V 
In this respect, model B did not succeed, and its projected 
distribution of HI in a trailing tail cannot be considered a 
satisfactory fit of the Magellanic Stream. 

Our models overestimate the integrated relative column 
densities of H I in the part of the Magellanic Stream located 
between the South Galactic pole and the far tip of the Stream. 
There is also no indication of the H I density decrease as we 
follow the Stream further from the Magellanic Clouds. In gen- 
eral, all of the models A, B, and C predict the trailing tail to be 
of higher H I column densities and extended well beyond the 
far tip of the Magellanic Stream. Such behavior is a common 
feature of pure tidal evolutionary models of the Magellanic 
System and it is a known drawback of omitting dissipative 
properties of the gaseous medium. 

Regarding the LSR radial velocity measurements along the 
Magellanic Stream by Briins et al. (2005 1, our models were 
able to reproduce some of their results. The far tip of the 
Magellanic Stream in model A reaches the extreme negative 
LSR radial velocity of -400 km s ' known from HI observa- 
tions. However, the highest negative LSR radial velocity does 
not drop below -350 km s ' for either prolate or nearly spher- 
ical halo configurations. Our previous discussion of various 
models of the Magellanic System denoted that the success- 
ful reproduction of the high-negative LSR radial velocity at 
the far tip of the Magellanic Stream is one of the most chal- 
lenging problems for such studies. Regarding our results, the 
importance of the coiTect LSR radial velocity profile along 
the Magellanic Stream was emphasized again. The absence 
of HI between LSR radial velocities of » -350 km s"' and 
~ -400 km s"' turned out to be the crucial factor decreasing 
the resulting fitness of examined evolutionary scenarios. Based 
on the study of debris of the Sagittarius dwarf galaxy. Law et 
al. ("2005") conclude that the velocity gradient along the trail- 
ing stream provides a sensitive evidence for the mass of the 
MW within 50 kpc and for the oblate shape of the galactic halo. 
However, the leading debris provide a contradictory evidence 
in favor of its prolate shape. 

5.4. Leading Arm 

Reproduction of the Leading Arm remains a difficult task for all 
the models of the Magellanic System that have been employed 
so far. Tidal models naturally place matter to the leading side 
of the System, towards the Galactic equator, as a result of the 
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tidal stripping also forming the trailing tail. However, neither 
the projected shape of the modeled leading structures nor the 
H I density distribution in the regions having an observational 
counterpart can be considered sufficient (see, e.g., Gardiner et 
al. [T994| . 

In every case A, B, and C, we were able to transport mat- 
ter to the area of the Leading Arm. Nevertheless, the pro- 
jected coverage of that region was more extended than what is 
observed. All the models contain a significant content of mat- 
ter spread from the leading edge of LMC across the Galactic 
equator, which has not been confirmed observationally. The 
model C reproduced the Leading Arm best. But the column 
density values of H I in C are overestimated and we also could 
not avoid an additional low-density envelope surrounding the 
structure (Fig.llSL 

6. Uncertainties in our modeling 

6.1. Missing pliysics 

To optimize performance of the GA, a computationally fast 
model of the Magellanic System is required. Therefore, 
complex N-body schemes involving self-consistent descrip- 
tions of gravity and hydrodynamics (see Bekki & Chiba 2005. 
Mastropietro et al. 2005 ) are discriminated. On the other hand, 
a correct description of physical processes dominating the evo- 
lution of the System remains a crucial constraint on the model. 

In Sect. 12. 21 we discussed the applicability of restricted N- 
body schemes on problems of galactic encounters and showed 
that they allow for modeling of extended streams and tails. 
Thus, we devised a restricted N-body code based on the nu- 
merical models by Murai & Fujimoto (1980 1 and Gardiner et 
al. (11994^ . The test-particle code interprets the observed large- 
scale structures such as the Magellanic Stream or the Leading 
Arm as products of tidal stripping in the Magellanic System. 

In addition to tidal schemes, ram pressure models have also 
been used in the previous studies on the Magellanic System 
(see Sect. 12.21 1. Compared to tidal schemes, hydrodynamical 
models allow for better reproduction of the H I column density 
in the Magellanic Stream, which is significantly overestimated 
if pure tidal stripping is assumed. Nevertheless, ram pressure 
models face difficulties concerning the clumpy structure of the 
Stream (see Briins et al. 2005i, which can hardly be reproduced 
by the process of continuous ram pressure stripping of the LMC 
gas as the Cloud moves through the halo of MW. On the other 
hand, the results of our study indicate that tidal stripping may 
produce inhomogeneous distribution of H I in the Magellanic 
Stream (see Fig.llO> due to the time variations of the tidal force 
exerted on the Clouds. Moreover, it is a significant drawback 
of ram pressure models that they constantly fail to reproduce 
the observed slope of the LSR radial velocity profile along the 
Magellanic Stream and especially the high-negative velocity 
tip of the Stream (Fig.|2}. Generally speaking, both classes of 
numerical schemes (tidal versus ram pressure) have introduced 
promising results into the research of the Magellanic System. 
Surprisingly, a study combining the effects of the tidal strip- 
ping due to the LMC-SMC-MW encounters with the influence 
of ram pressure stripping is still missing. The recent papers ei- 



ther do not focus on the large scale structure of the Magellanic 
System (e.g., Bekki & Chiba 2005 ), or ommit the role of grav- 
itational encounters by neglecting the LMC-SMC interaction 
and restricting the studied evolutionary period of the System 
(see, e.g., Mastropietro et al. 2005)- 

Both families of models suffer from serious difficulties 
when modeling the Leading Arm. The tidal models are able 
to place matter into the area of the Leading Arm, since the cre- 
ation of a trailing tail (the Magellanic Stream) is naturally ac- 
companied by the evolution of a leading stream (for details see 
Toomre & Toomre 1972 ). However, neither the projected shape 
of the modeled Leading Arm nor the local column density of 
HI satisfy the observations. Implementation of a hydrodynam- 
ical scheme may improve the results concerning the spatial ex- 
tent of the tidal Leading Arm and its H I column density profile, 
due to ram pressure exerted by the MW gas on the leading side 
of the Clouds. Thus, it is another argument for further consid- 
eration to be given to a hybrid tidal-i-ram pressure model of the 
LMC-SMC-MW interaction. 

From a technical point of view, employing even a simple 
formula for ram pressure stripping would introduce other pa- 
rameters including structural parameters of the distribution of 
gas in the MW halo and a description of the gaseous clouds in 
the LMC and SMC. It would increase the dimension of the pa- 
rameter space of the interaction and complicate the entire GA 
optimization process. 

6.2. Mass and shape evolution of the Magellanic 
Clouds 

The dark matter halo of the MW is considered axisymmetric 
and generally flattened in our model. It is a significant improve- 
ment over previous studies of the Magellanic System that as- 
sumed spherical halos only. We were able to investigate the in- 
fluence of the potential flattening parameter q on the evolution 
of the Magellanic System. However, both the mass and shape 
of the MW DM halo were fixed for the entire evolutionary pe- 
riod of 4 Gyr. 

We did not take into account possible changes in mass and 
shape of the Clouds. Shape modification might become impor- 
tant for very close LMC-SMC encounters that are typical for 
the models with nearly-spherical MW DM halos. Penarrubia et 
al. J2004> demonstrated that a relative mass-loss of a satellite 
galaxy moving through an extended halo strongly differs for 
various combination of its orbital parameters, shape, and the 
mass distribution of the halo, and cannot be described reliably 
by a simple analytic formula. 

7. Summary and conclusions 

We performed an extended analysis of the parameter space for 
the interaction of the Magellanic System with the Milky Way. 
The varied parameters cover the phase space parameters, the 
masses, the structure, and the orientation of both Magellanic 
Clouds as well as the flattening of the dark matter halo of the 
MW. The analysis was done by a specially adopted optimiza- 
tion code searching for a best match between numerical models 
and the detailed HI map of the Magellanic System by Briins et 
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al. <2005> . The applied search algorithm is a genetic algorithm 
combined with a code based on the fast, but approximative re- 
stricted N-body method. By this, we were able to analyze more 
than 10^ models, which makes this study one of the most ex- 
tended ones for the Magellanic System. 

In this work we focused especially on the flattening of the 
MW dark matter halo potential q. The range 0.74 < q < 1.20 
was studied. It is equivalent to the interval of the density flat- 
tening 0.31 <qp< 1.37 (see Eq.|3}. 

We showed that the creation of a trailing tail (Magellanic 
Stream) and a leading stream (Leading Arm) is quite a common 
feature of the LMC-SMC-MW interaction, and such structures 
were modeled across the entire range of halo flattening values. 
However, important differences exist between the models, con- 
cerning density distribution and kinematics of H 1, and also dy- 
namical evolution of the Magellanic System over the last 4 Gyr. 
In contrast to Murai&Fujimoto ( 1980l, Gardiner et al. ( 1994 1, 
or Lin et al. {1995 1, the Clouds do not have to be gravitationally 
bound to each other for the entire evolutionary period to pro- 
duce the matter distribution that is in agreement with currently 
available H 1 data on the Magellanic System. 

Overall agreement between the modeled and observed dis- 
tribution of neutral hydrogen in the System is quantified by 
the fitness of the models. The fitness value is returned by a 
FF, that performs a very detailed evaluation of every model 
(AppendixiBl. Analysis of fitness as a function of the halo flat- 
tening parameter q indicates that the models assuming oblate 
DM halo of MW (model A) allow for better satisfaction of 
HI observations than models with other halo configurations. 
Analysis of Fig.|4]does not indicate a drop in the value of fit- 
ness as the flattening parameter q decreases. It suggests that 
models of a quality comparable to the fitness of model A may 
exist for even more oblate configurations of the MW DM halo. 
Unfortunately, the logarithmic potential defined by Eq.^ can- 
not be used for q < 1/ y/2 due to negative density on the axis of 
symmetry. Thus, the result represented by Fig.l^established a 
strong motivation for further extension of the parameter study 
using a more general and realistic model of the Galactic halo 
(e.g., a tri-axial distribution of matter or the flattening parame- 
ter depending on position). Nevertheless, the above mentioned 
problem does not alter our result preferring flattened halo con- 
figurations and discriminating nearly-spherical shapes within 
the galactocentric radius of « 200 kpc containing typical orbits 
of the Clouds and so analyzed in this study. 

We did not involve surveys of stellar populations in the 
Magellanic System in the process of fitness calculations. This 
is due to the nature of test-particle models that do not allow for 
the distinction between stellar and gaseous content of studied 
systems. However, we still have to face one of the most inter- 
esting observational facts connected to the Magellanic Clouds 
- the absence of stars in the Magellanic Stream (Van den 
Bergh l200()t - because both the LMC and SMC contain stellar 
populations, and so every structure emanating from the Clouds 
should be contaminated by stars. It is an additional constraint 
on the models. It cannot be involved in fitness calculation be- 
cause of the limits of our numerical code, but has to be taken 
into account. 



Stellar populations of the SMC are very young and the mass 
fraction in the form of stars is extremely low. Our models show 
that the evolution of the Magellanic Stream has been lasting 
2 Gyr at least (model B). Thus, the fraction of matter in the 
Magellanic Stream, that is of SMC origin, was torn off before 
significant star formation bursts occurred in the SMC, and stars 
should not be expected in the Stream. Nevertheless, we found 
both LMC and SMC matter in the Magellanic Stream for ev- 
ery model of the System. Similarly to the case of the SMC, if 
the LMC star formation activity was increased after the matter 
transport into the Magellanic Stream was triggered, stars would 
naturally be missing in the Stream. Such a scenario is doubtful 
however Observational studies argue for a massive star forma- 
tion burst started in LMC at ^ -3 Gyr (Van den Bergh 20^, 
which is rather close to the age of the Magellanic Stream, as 
indicated by our models (see Sect. 15. 3> . Our results concerning 
the LMC and SMC orbits introduced an acceptable solution 
to the problem of missing stars. We showed that the evolution 
of the Clouds in aspherical MW DM halos (models A and C) 
does not lead to extremely close encounters disturbing inner 
parts of the LMC disk (rdisk < 5 kpc). Since the distribution 
of gaseous matter in galaxies is typically more extended than 
the stellar content, the Magellanic Stream matter coming from 
outer regions of the Clouds does not necessarily have to contain 
a stellar fraction. 

Previous discussion of stellar content of the Magellanic 
System supports discrimination of the configurations with 
nearly spherical halos (model B) that was discovered by the 
GA search. On the other hand, many papers on the dynami- 
cal evolution of the Magellanic Clouds dealing with a spher- 
ical MW halo (Murai&Fujimoto 1980 Gardiner et al. [T^ 
Bekki&Chiba 2005) argue that the observed massive LMC 
star formation bursts 3 Gyr ago was caused by close LMC- 
SMC encounters. Our model B shows close approaches of the 
Clouds Ar a: lOkpc at around the mentioned time. For as- 
pherical halos, such encounters do not induce the formation 
of particle streams. However, close LMC-MW and SMC-MW 
encounters appeared to be efficient enough to trigger massive 
matter redistribution in the System leading to formation of the 
observed structures. Then, they could also be responsible for 
the triggering of star bursts. 
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Appendix A: Dynamical friction 

If the distribution function in velocity space is axisymmetric, 
the zeroth order specific friction force is (Binnev ri977> : 



It was shown by Gardiner et al. J1994t and 



F' 



2 ^pi^iR, z)G^Ms y/l-ey^\n A 
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2 V2^pl(/?, z)G^Ms ^Jl - e,,^ In A 
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(A.l) 



(A.2) 



where i - x,y and (ctr, ct,) is the velocity dispersion ellipsoid 
with ellipticity = 1 - (cr,/crR)^, InA is the Coulomb log- 
arithm (Chandrasekhar il943< of the halo, Ms is the satellite 
mass, and 
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where (vr, w-) are the components of the satellite velocity in 
cylindrical coordinates. 

Appendix B: Fitness function 

The behavior of the 3 D test-particle model of the Magellanic 
System is determined by a large set of initial conditions and 
parameters that can be viewed as a point (individual) in the 
system's high-dimensional parameter space. In the case of our 
task, the fitness of an individual means the ability of the nu- 
merical model to reproduce the observed H I distribution in the 
Magellanic Clouds if the individual serves as the input param- 
eter set for the model. It is well known that proper choice for 
FF is critical for the efficiency of GA and its convergence rate 
to quality solutions. After extended testing, we devised a three- 
component FF scheme. To discover possible unwanted depen- 
dence of our GA on the specific choice for the FF, both of the 
following FF definitions were employed: 



FFa — FFi ■ FF2 ■ FF3, 



FFb = 



1=3 

2 c,- ■ FFi 

i=i 

(=3 

!=1 



(B.l) 



(B.2) 



where the components FFi, FF2, and FF3 reflect significant 
features of the observational data and ci = 1.0, C2 = 4.0, and 
C3 = 4.0 are weight factors. Both the FFs return values from 
the interval (0.0,1.0). 

FF compares observational data with its models. To do that, 
resulting particle distribution has to be treated as neutral hydro- 
gen and converted into H I emission maps for the defined radial 
velocity channels. In the following paragraphs we briefly intro- 
duce both the observed and modeled data processing. 



Gardiner & Noguchi J1996t that the overall HI distribu- 
tion in the Magellanic System (Magellanic Stream, Leading 
Arm) can be considered a tidal feature. Following that result, 
an elaborate scheme of the original data (Briins et al. I2005> 
manipulation was devised to emphasize large-scale features 
of the Magellanic HI distribution on the one hand, and to 
suppress small-scale structures on the other hand, since 
they originate in physical processes missing in our simple 
test-particle model. The observational data are stored in the 
Flexible Image Transport System (FITS) format, which lets us 



Magellanic Stream 




Fig. B.l. The figure depicts the original 3D HI data cube 
by Briins et al. f2005'l (upper plot) together with the re- 
sulting data after median and Fourier filtering. Both im- 
ages offer 3D visualization of the column density isosurface 
Shi =0.2- 10'^— -2 



cm 



apply standard image processing methods naturally. A Fourier 
filter was selected for our task. It represents a frequency 
domain filter, and so it allows for an excellent control over 
the scale range of the image's structures to be conserved or 
filtered out. We removed the wavelengths below the limit of 
« 10° projected on the sky-plane. The performance of Fourier 
filters suffers from the presence of abrupt changes of intensity, 
such as edges and isolated pixels. To enhance the efficiency 
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of frequency filtering, it was preceded by an application of 
a spatial median filter to smear the original image on small 
scales. Subsequently, the HI column density is normalized. 
The resulting 3D HI column density data cube together 
with the original data by Briins et al. (2005) can be seen in 
Fig. IB. II To compare the modeled particle distribution with 
H I observations, we convert the distribution to a 3 D FITS 
image of column densities that are proportional to particle 
counts, since all the test-particles have the same weight factor 
assigned. Then, we have to interpolate missing data, which is 
due to a limited number of particles in our simulations. Finally, 
the column density is normalized to the maximal value. 

After discussing the data processing and manipulation, we 
will introduce the individual FF components FFi, FF2, and 
FF3. 

B.1. FFi 

The observed HI LSR radial velocity profile measured 
along the Magellanic Stream is a notable feature of the 
Magellanic System. It shows a linear dependence of LSR 
radial velocity on Magellanic Longitude, and a high neg- 
ative velocity of -400 km s^' is reached at the Magellanic 
Stream's far tip (Briins et al. '2005 1. From the studies by 
Murai & Fujimoto ( 1980j and Gardiner et al. ( 1994i, and from 
our modeling of various Magellanic evolutionary scenarios, 
we know that the linearity of the Magellanic Stream veloc- 
ity profile shows low sensitivity to the variation of the initial 
conditions of the models. On the other hand, the slope of the 
LSR radial velocity function is a very specific feature, espe- 
cially strongly dependent on the features of the orbital motion 
of the Clouds. Therefore, it turned out to be an efficient ap- 
proach to test whether our modeled particle distribution was 
able to reproduce the high negative LSR redial velocity tip of 
the Magelanic Stream. Then, the first FF component FF\ was 
defined as follows: 



FF, = 



1 



1 -I- 



(B.3) 



where v°'°^ and v™'^ are the minima of the observed LSR ra- 

min min 

dial velocity profile of the Magellanic Stream and its model, 
respectively. 



The second FF component analyzes whether there is a mod- 
eled H I emission present at the positions and LSR radial veloc- 
ities where it is observed. Thus, we measure the relative spatial 
coverage of the System observed in H I emission by the mod- 
eled matter distribution for every LSR radial velocity channel. 
No attention is paid to specific H I column density values here. 
We only test whether both modeled and observed emission is 
present at the same pixel of the position-velocity space. It can 
be expressed as 
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where pix°j^ e {0, 1 ) and pixfj^'^ e {0, 1 ) indicate whether there 
is matter detected at the position [i, j, k] of the 3 D data on the 
observed and modeled Magellanic System, respectively. - 
64 is the number of separate LSR radial velocity channels in 
our data. {N^ ■ Ny) - (64 • 128) is the total number of positions 
on the sky-plane for which observed and modeled H I column 
density values are available. 

This binary comparison between the observed and modeled 
data introduces a problem of pure noise pixels present in the 
observed data cube because they posses the same weight as the 
other data, despite their typically very low intensity. However, 
our treatment of the original high-resolution data by Briins et 
al. (2005 1 involves spatial median filtering. It smears abrupt in- 
tensity changes and removes isolated pixels, which handles the 
problem of pure-noise data pixels naturally. The subsequent 
Fourier filtering decreases the data resolution significantly, and 
that also strongly suppresses the influence of original noise pix- 
els. 

As the last step we compare the modeled matter density 
distribution to the observation. To do that, both modeled and 
observed H I column density values are scaled relative to their 
maxima to introduce dimensionless quantities. Then, we get 
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where cr°^^, cr™!"^ are normalized column densities measured at 

i]k ' t jk 

the position [j, k] of the i-th velocity channel of the observed 
and modeled data, respectively. 



B.2. FF2 and FF3 

The FF components FF2 and FFj compare the observed and 
modeled HI column density distributions in the Magellanic 
System for 64 separate LSR radial velocity channels of width 
Av = 13.2 km s '. For every velocity channel, H I column den- 
sity values are available for (64 ■ 128) pixels covering the entire 
System. The above introduced 3 D data was obtained by modi- 
fication of the original high-resolution H I data-cube by Briins 
et al. ( I2005> . Since the test-particle model is not capable of re- 
producing small-scale features of the explored system, filtering 
and reduction of resolution of the original data were necessary 
prior to its use for the purpose of our GA search. 



